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Abstract 

We discuss the relevance of studying ecology within the framework of Com- 
plexity Science from a statistical mechanics approach. Ecology is concerned with 
understanding how systems level properties emerge out of the multitude of interac- 
tions amongst large numbers of components, leading to ecosystems that possess the 
prototypical characteristics of complex systems. We argue that statistical mechan- 
ics is at present the best methodology available to obtain a quantitative description 
of complex systems, and that ecology is in urgent need of "integrative" approaches 
that are quantitative and non-stationary. We describe examples where combining 
statistical mechanics and ecology has led to improved ecological modelling and, at 
the same time, broadened the scope of statistical mechanics. 
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1 Introduction 

Does the concept of "Complexity" bear any specific meaning or is it just synonymous 
with complicated and yet not comprehended phenomena? We will argue that it is pos- 
sible and useful to use the term "Complexity Science" in a specific and reasonably well 
denned way. It is useful because a number of common trends and implications become 
clear when a phenomenon is classified as part of "Complexity Science". The science of 
complexity emphasises the interactions between components. It stresses that compo- 
nents, most often, are heterogeneous and evolve in time. Complexity is concerned with 
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the emergent properties at systems level originating from the underlying multitude of 
microscopic interactions. 

In an attempt to make our discussion more clear we will immediately describe the way 
we use some terms central to our exposition. We hurry to stress that these descriptions 
are not meant to be exhaustive final philosophical definitions, but rather intended to 
lower the risk of misunderstanding when we deal with terms frequently used to mean 
different things by different people. And now our specifications. 

Complex Systems consist of a large number of interacting components. The 
interactions give rise to emergent hierarchical structures. The components of the 
system and properties at systems level typically change with time. A complex 
system is inherently open and its boundaries often a matter of convention. 

Statistical Mechanics seeks to understand how properties at systems level emerge 
from the level of the system-components and their interactions. This often involves 
the application of probability theory, and a number of mathematical techniques. 
Throughout, we draw a distinction between statistical mechanics and statistical 
physics. The latter is mainly concerned with the microscopic foundation of ther- 
modynamics and, e.g., phenomena such as phase transitions and superconductivity. 
We consider here statistical mechanics as a mathematical methodology, which can 
be applied to many different sciences including economics, population biology and 
sociology, to name a few. 

Statistical mechanics is a powerful transdisciplinary methodology for the study of 
emergent phenomena at a macroscopic level caused by the many interactions taking place 
at a microscopic level. It provides a framework within which it is possible to encapsulate 
the myriad of degrees of freedom of a system at a microscopic level, into just a few degrees 
of freedom at a macroscopic level. In its current form statistical mechanics does not hold 
all the answers for all the complex systems, however, we argue that it is at present the 
best methodology available to obtain a quantitative description of complex systems. By 
systematically applying it to fields outside its traditional range of application in physics, 
statistical mechanics can be developed further, in addition to simultaneously contributing 
to the understanding of those fields, such as ecology. The importance of this feedback loop 
cannot be overestimated. It can also provide a starting point for the possible development 
of new mathematical techniques. 

Along these lines, the research programme in search for the "laws" of ecosystems 
described by J0rgensen and collaborators [T|, looks into finding a rigorous set of laws 
that govern the dynamics at the macro-level. This is a first attempt into establishing a 
methodology for ecological complexity. At the moment the analysis is mainly qualitative, 
and we suggest that the second step towards that goal would be the implementation 
of techniques from statistical mechanics, in order to obtain a rigorous mathematical 
formalism and modelling. 

The paper is organised as follows. In the next section, we will for concreteness il- 
lustrate our arguments by briefly describing a complexity inspired model of evolutionary 
ecology called the Tangled Nature model. This will allow us to demonstrate how macro- 
evolution can be modelled as emerging from the interacting micro-evolution consisting 
of individual organisms influencing each other and undergoing reproduction which is 
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prone to mutation. We will discuss feedback, emergence, network structures, and the 
intermittent temporal mode of macro-evolution in contrast to the steady smooth pace 
of dynamics at the level of individuals. We will also briefly touch on the modelling of 
ecological observables such as the Species Abundance Distribution, Species Area Laws 
and the relationship between interaction and diversity. 

In the Discussion Section we will mention two examples of research presented at the 
Symposium on Complexity, Collective Effects and Modelling of Ecosystems: formation, 
function and stability at the Beijing Eco Summit 2007. These examples illustrate how 
a complexity science viewpoint may shape the approach of ecological research projects. 
The first example is John Crawford's contribution on "The Self-organisation of life in 
Earth" . This work looks at the soil-microbe system, and at the development of models on 
evolutionary ecology, that can be applied to this dynamical system. The second example 
is Cedric Gaucherel's work on "Theoretical analysis of dynamic patchy landscapes", which 
looks at landscape models constructed within the framework of statistical mechanics. 

2 Tangled Nature Model 
2.1 Description of model 

The Tangled Nature model is defined at the level of interacting individuals. It is an 
attempt to identify possible simple mechanisms behind the myriad of complicated inter- 
actions, feedback loops, contingencies, etc., as one moves from the short time reproductive 
dynamics at the level of individuals, to the long time systems level behaviour. The strat- 
egy is to keep the model sufficiently simple to enable analysis and to pinpoint the details 
or assumptions in the model that are responsible for the specific behaviour at the systems 
level. One major concern of the model has been to understand how the smooth continu- 
ous pace of the reproductive dynamics at the level of individuals, can lead to intermittent 
or punctuated dynamics at the level of high taxonomic structures. To be able to address 
such issues, the model considers individuals as represented by a single sequence with 
individual number a, denoted by S a = (Sf , S% , belonging to a sequence space 

iS, where all £■* = ±1. These sequences undergo simple reproduction during which a 
given sequence duplicates itself, and while this happens, components of the sequence 
may mutate, represented by the off-spring having a different sign from the mother, i.e. 
5*7 = — S" , where 7 denotes the index for the daughter, and a the one for the mother. 
The aim of the model is to understand the macro-dynamics emerging at the systems 
level. This is done by analysing the dynamics of the occupancy in this sequence or type 
space. Taxonomic structures, such as species formation, emerge as aggregations in the 
density of occupied sites n(S,t) in the type space. This is very much in accordance with 
Mallet's definition of species [2J. A species will be identified as a local peak in the density 
n(S,t), and species formation will correspond, e.g., to the splitting of such a peak into 
two peaks. Macroscopic ecological measures such as Species Abundance Distributions are 
derived from the structure of n(S,t). 

Let us now sketch the mathematical details of the model. For in depth studies of the 
model, please refer to [SI IH El E]- The size of the type space is set by the length, L, 
of the sequences; a typical value used is L = 20 leading to about one million different 
genotypes. The sites in the genome space are supposed to represent all possible ways 
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of constructing a 'genome'. Many sequences may not correspond to viable organisms. 
The viability of a genotype is determined by the evolutionary dynamics. All possible 
sequences are available for evolution to select from. We will see that a natural species 
concept arises from the dynamics, in which each species is separated in genotype space. 

The system consists of N(t) individuals, and a time step consists of one annihilation 
attempt followed by one reproduction attempt. A reproduction event is successful with 
varying probability p ff, defined below, and an annihilation attempt is successful with 
constant probability pun- The killing probability is considered a constant independent 
of type for simplicity. It would obviously be more realistic to let pkm depend on the 
type of the individual considered. However, this does not change the overall behaviour at 
systems level. One generation consists of N(t)/pkm time steps, which is the time taken 
(on average) to kill all currently living individuals. The dynamics leads to a population 
size which remains nearly constant on short timescales. The individuality of the specific 
types, or sequences, is given by their ability to reproduce. Since we are interested in 
the collective, or complexity, aspects of evolution, the Tangled Nature model stresses the 
mutual influence amongst different types of organisms. This is done by assuming that 
each individual of type S is able to reproduce, when selected for reproduction, with a 
probability p //(S, t) that depends on the sequence S and the configuration of other types 
in the type space. The reproduction probability, p ff, is determined by a weight function 
H{S a ,t): 

^ (sa ' t] = W) (jy j(s °' sms ' t] ) ~ (i) 

where c controls the strength of the interaction (large c means a large interaction), N(t) 
is the total number of individuals at time t, the sum is over the 2 L locations in S and 
n(S,t) is the number of individuals (or occupancy) at position S. Two positions S a and 
S b in genome space are coupled with fixed but random strength J(S a , S b ) which can be 
either positive, negative or zero. This link exists (in both directions) with probability 9, 
i.e. 9 is simply the probability that any two sites are interacting. If the link exists, then 
J(S a ,S b ) and J(S 6 , S a ) are both generated randomly and independently, and such that 
they belong to ( — 1, 1). To study the effects of interactions between species, we exclude 
self-interaction so that J(S a , S a ) = 0. 

The conditions of the physical environment are simplistically described by the term 
fiN(t) in Eq. ([T]), where /i determines the average sustainable total population size, i.e. 
the carrying capacity of the environment. This is an example of how the question of 
the openness and "surroundings" of ecosystems arises in a natural way in the present 
statistical mechanics like formalism. An increase in /i corresponds to harsher physical 
conditions. We use asexual reproduction consisting of one individual being replaced 
by two copies mimicking the process of binary fission seen in bacteria. We allow for 
mutations in the following way: with probability p mu t per gene we perform a change of 
sign S" — > — S'f during reproduction. Successful reproduction occurs with a probability 
per unit time, p ff(S a ,t) e [0, 1], given by 

„ fS q f) exp[ff(SV)] 

Poff{S '^-l + exp^SV)]- (2) 

This function is chosen for convenience, since the specific functional form has no effect on 
the dynamics of the model - any smoothly increasing function that maps H(S a , t) to the 
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interval [0, 1] will do. Let us mention that this basic quantity is deliberately taken by the 
Tangled Nature model to be a context dependent reproduction probability rather than a 
fitness function. One reason why this is done is to try to avoid some of the dangers and 
subtleties inherent to the fitness concept [7]. 

Eq. ([T]) can be understood as the interaction of an individual with all the others, 
with a term fiN which determines the total population and controls fluctuations. The 
interaction strength c gives the magnitude of the total interaction. We can tune the 
effective 'resource' density (and hence the population density) with the parameters \i. The 
total population remains approximately constant over ecological timescales (and actually 
increases over evolutionary timescales). Setting self-interaction to zero is equivalent to 
considering that all types interact equally with their own species (one can rescale pkm 
and n to accommodate this). This constraint is imposed in order to focus on the effect 
of interactions between different types. To study the relation between diversity and the 
strength of the interactions, there is a version of the model in which different strengths 
of intra-specific interactions are included [Hj. Obviously it is impossible to design the 
details of the interaction matrix J(S a , S b ) in a realistic way. What can be accomplished 
is to study qualitative questions, such as what is the effect of having very few interaction 
links between sequences compared with many interaction links Or one can address 
the effect of correlations in the allowed interactions [9]? 

After a short transient period the initial state becomes irrelevant. There are two 
very different initial conditions that consist in placing the entire population at time zero 
in: i) one position in type space, or ii) on random positions, i.e. a random collection 
of initial types. Since both configurations are badly adapted to the interaction matrix 
J(S a , S b ), in both cases the population will typically collapse to one single position in type 
space. Eventually the population size will have decreased enough to make the —p,N(t) 
term sufficiently small to allow p Q ff to grow to a value that ensures a non- vanishing 
reproduction rate. When this happens the population will, due to mutations, start to 
spread out from its initial position into the surrounding genotype space. And as this 
happens, natural selection will ensure that only certain configurations of occupied sites 
are viable. These are configurations for which the mutual interactions between the types 
lead to off-spring probabilities that, for a significant part of the occupied types, are able 
to balance the killing probabilities, i.e. p a ff(S,t) = pkm for some set of types S. 

The dynamics in type space is characterised by a two-phase switching, consisting 
of long periods of relatively stable configurations (quasi-Evolutionary Stable Strategies 
or q-ESSs) (Fig. [T]) interrupted by brief spells of reorganisation of occupancy called 
transitions. Transition periods are terminated when a new q-ESS is found, as discussed 
in [3J. The intermittent macro-dynamics is not in a stationary state. When one considers 
very many realisations of the dynamics it turns out that the transition rate between 
q-ESS decreases with the age of the system [6]. This happens because selection is able 
to pick out configurations in type space that tend to possess more beneficial links (i.e. 
positive J(S a , S b ) bonds) than is the case between a randomly selected set of types. We 
consider this directedness of the long time systems level dynamics to prototypical of 
complex systems [TO] . 
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Figure 1: Intermittent evolution of the occupancy in type space. Time, measured in 
generations, is along the x-axis. The about 10 6 different types are labelled up along the 
y-axis. Whenever a type is occupied a dot is placed at its label. Long stretches of parallel 
lines indicate epochs during which the main composition in type space remains essentially 
the same. Figure courtesy of Matt Hall. 

2.2 Emergent Ecological Measures 

As we move from the level of individuals to the systems level, 'species' can be well 
defined as the highly occupied genotype points called 'wildtypes', which are separated in 
genotype space. Each wildtype is surrounded by a 'cloud' of mutant genotypes with low 
occupancy. Thus we can take a natural definition of diversity: the number of wildtypes 
in the system. It is interesting to study how the Species Abundance Distribution (SAD) 
depends on the assumed properties of the interaction matrix J(S a , S b ). It was found by 
|4J that the often observed log-normal shape of the SAD is reproduced by the evolutionary 
dynamics of the Tangled Nature model under certain conditions. Namely, when each type 
is potentially able to interact with a large number of other types. In this case the adapted 
configurations consist of populations of species that form one large interconnected cluster 
and the SAD evolves with time towards a log-normal like form. If J(S a , S b ) only allows 
a type to interact with few other types (i.e. few non-zero elements in the J matrix), the 
population in the type space splits up into separate groups and the SAD doesn't develop 
a form resembling a log-normal distribution. Considered from this perspective the SAD 
might be thought of as containing information about the properties of the network of all 
possible interactions between organisms. 

Let us focus on the properties of the network of interactions^ of extant species. The 
evolutionary dynamics performs a collective adaptation on the co-existing types in type 
space. As selection and adaptation act generation after generation, a subset of sites 
in type space becomes occupied. This subset is selected such that the mutual interac- 
tions allow each of the extant species to counter balance the depletion of its population, 
caused by death (pkui) and mutations (p m ut), by a sufficiently large off-spring production 
(p a ff(S,t)). The network of interactions between these co-existing types possesses some 
interesting emergent properties. The typical coupling strength between extant types is 
more mutualistic than the coupling between arbitrary types S a and S b , chosen at random 
in type space irrespectively of the types being extant or not [U [9]. This effect is signif- 

*The nodes of the network under consideration, consist of occupied positions in type space. There is 
an edge between two nodes S a and S h if the two types interact, i.e. if J(S a , S b ) or J(S b , S a ) is non-zero. 
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icantly bigger when the coupling matrix J(S a ,S b ) is correlated for sites S a and S 6 that 
reside in the same vicinity of type space. A correlated coupling matrix is more realistic, 
since it corresponds to assuming that similar organisms have a certain similarity in the 
way they interact with the surrounding ecosystem. The degree distribution of the net- 
work of extant species is sensitive to the amount of correlations imposed on type space. 
When correlations are present, we typically observe exponential degree distributions of 
the network of interactions between extant types. In contrast, un correlated interaction 
matrices lead to binomial degree distributions, as it is observed in networks where edges 
are placed at random. This example indicates how some properties at systems level may 
be caused by generic mechanisms for emergent collective behaviour. 

A similar situation is encountered when the qualitative behaviour of spatial properties 
is investigated. Spatial aspects are obviously of the greatest ecological importance. A 
simple quantity to start out with is the Species Area Relation (SAR). By placing a copy 
of the Tangled Nature model on each site of a two dimensional lattice, one can make 
a simplistic model combining evolutionary dynamics with spatial dispersion. Such a 
model was studied in [TT] and a power-law SAR is observed. The evolutionary dynamics 
produces a high degree of spatial diversity even when the same type space is placed on 
each site of the spatial lattice. 

Complexity science stresses that the interaction between the components is responsible 
for the emergent properties at systems level. Sometimes for tractability reasons, models 
might oversimplify the components compared with reality, and still it happens that such 
simple models are able to capture certain qualitative aspects. An attempt in this direction 
was made in [8], where the relationship between interactions amongst different types 
and the diversity of types was discussed. The inspiration behind this study came from 
molecular evolution experiments on E. coli, in which the relation between fitness plasticity 
and diversification was addressed [T2]. The model used a version of the Tangled Nature 
model in which each type was assigned an amount of self-interaction. This was achieved by 
supplementing the weight function in Eq. (CQ) by an extra term proportional to an intrinsic 
fitness E(S a ). It was found that diversity rapidly increased when the typical interaction 
strength, set by the parameter c in Eq. (JTJ), exceeded a certain value, determined by the 
properties of the intrinsic fitness. 

3 Discussion 

We have proposed that complexity science can be seen as a coordinated attempt to un- 
derstand how emergent collective behaviour at systems level arises due to the multitude 
of interactions between the components. From this perspective we have argued that 
complexity science offers a particularly relevant approach to ecology. Above we tried 
to illustrate our point of view with some theoretical examples taken from the Tangled 
Nature model's study of evolutionary ecology. Let us now conclude with a couple of 
examples that are closely related to observations, where we believe that the complexity 
science's perspective has made a difference. 

Our first example is the soil-microbe system. In the study of such systems, there is an 
urgent need to develop models on evolutionary ecology that integrate function and diver- 
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sity, and that are dynamical p3] , at present many of the current models for soil are static. 
The soil-microbe system is an extremely rich and intricate system that has not yet been 
fully understood, and moreover is of great relevance to agriculture, waste management 
and the water industry to mention a few [14J. Crawford and collaborators believe that 
any model describing the system should integrate biochemistry and biophysics, since from 
the interactions between the biotic and abiotic factors, the structure, the functionality 
and the dynamical behaviour of the soil emerge. 

The soil- microbe complex is a self-organised system capable to adapt, therefore, within 
the many different approaches that currently exist [13], they suggest the development of 
models of evolutionary ecology that have the same aspects as those described by the 
Tangled Nature model: evolutionary processes, population dynamics, feedback loops, in- 
teractions, etc. In addition, the ecosystem functioning needs to be included, since it is the 
relation between diversity and function that is mostly ignored in other soil models. Ex- 
perimental techniques to measure this connection need therefore to be developed. They 
argue that the ecosystem-level behaviour of the soil-microbe system is the outcome of the 
behaviour at the organism's level, which is natural if the system is classified as a complex 
system. In addition, this viewpoint stresses the importance to include evolutionary pro- 
cesses when looking at ecology, since these are crucial for the understanding of ecological 
function. For example, in the soil-microbe system, the activity of the microbes changes 
the structure of the soil by affecting its rates of oxygen diffusion and porosity, while the 
substrates in the soil affect the activity of the microbes. 

Crawford and collaborators have identified important properties of soil systems by 
applying statistical mechanics to their research [15J, however, they also recognise the 
need of an extended methodology, where interdisciplinarity is crucial. They emphasise 
that there is no unique discipline that is able by itself to understand the soil [Tj] . 




Our second example is the modelling of landscapes. Models explaining and being able 
to predict the shape of landscapes are extremely important to prove ecological hypothe- 
ses, and for the implementation and development of market and land planning policies. 
Gaucherel and collaborators [TB] . argue that in many systems, the most relevant factors 
causing the dynamical changes of landscapes are human driven. For this reason they 
urge for the development of models that integrate biophysical and socio-economical fac- 
tors. They propose a generic modelling platform: "LI", that can be used to look at 
the patterns resulting from specific processes. This can therefore be applied as a tool to 
assess environmental policies and technological implementations at different landscape 
scales. The platform simulates the dynamical evolution of a landscape as a result of the 
feedbacks and interactions between the elements composing the landscape. Following the 
methodology of statistical mechanics, the objects are modelled by introducing only the 
relevant aspects that give rise to such emergent structures, and not by parametrising all 
of their degrees of freedom. 

In their approach, Gaucherel and collaborators stress the crucial role of feedback and 
scales in ecological systems. On the one hand, the landscape is an emergent structure, 
product of the interactions between the many different components, see for example the 
role of farm systems in human driven landscapes [IE]. On the other, their models give 
rise to hierarchical structures that feed back into the system. For example, the landscape 
itself determines important aspects of habitats and ecosystems at different scales. This is 
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illustrated in [T7], where the authors look at the relationship between the characteristics 
of the habitat given by the landscape, and particular characteristics of the inhabitant 
species, such as their spatial distribution, their morphology, etc. Using their framework, 
one can investigate at different locations, times and scales, the ecological relationships. 
They apply their model to look at the link between agricultural activities, landscape 
shape and some characteristics of carabid beetles, such as their abundance and their 
body size. They proved that some correlations are only valid for certain specific locations 
and scales, contrary to what was believed using other techniques unable to give spatial 
and scaling precision. 

An explicit outline of how the methodology from statistical mechanics is implemented 
in this approach to landscape modelling can be found in [TH] . There the authors construct 
a neutral model for patchy landscapes using the Gibbs process to describe the interac- 
tions between the different components of the landscape. A neutral model is a model that 
simulates the landscape properties and patterns that are not the outcome of a particular 
ecological process. Therefore, using these models a distinction can be drawn between 
structures caused by random processes and those obtained through real processes. In 
addition, such models give rise to virtual landscapes that can be used to study real mo- 
saics, such as forest landscapes [H]. For this or other specific applications, the neutral 
model is implemented in the LI platform mentioned above, where the particularities of 
the system can be entered. 

The project to establish a systems perspective on ecology as laid out in the book A New 
Ecology by J0rgensen and collaborators in pQ can, in our opinion, be seen as a prototypical 
example of the objectives aimed at when taking inspiration from the methodology of 
statistical mechanics and applying it to complex systems. J0rgensen et al. argue that 
laws at the emergent systems level may exist and the authors suggest a list of laws 
they believe ecosystems obey. They make clear that it is a grand task to identify these 
laws and that their list is to be thought of as a starting point. If we embed these laws 
within the statistical mechanics framework, they describe the expected properties of a 
complex system. For example, J0gensen et al. mention that ecosystems have openness, 
connectivity, complex dynamics and that their dynamics is directed. These properties 
are totally in agreement with the features of complex systems highlighted by applying 
methods from statistical mechanics to their analysis. 

In addition, within the Tangled Nature model, interconnectedness and interaction be- 
come a focal point of the description when one thinks in terms of the emergent networks 
of interactions between extant species. Furthermore, within the mathematical formula- 
tion of dynamical systems, there is a term encoding the carrying capacity and resources, 
which represents an open system. Ecosystems are therefore correctly taken into account 
as open within this formalism. The Tangled Nature model was formulated by including 
what appears to be minimal assumptions for the dynamics; namely reproduction prone 
to mutations and livelihoods of the individual types that are influenced by other coex- 
isting types. As a result the model produces a slowly adapting non-stationary directed 
dynamics at the macroscopic systems level. This is certainly in agreement with prop- 
erties of ecosystems encoded in the form of a law by J0rgensen et al. and, moreover, 
it appears to be in agreement with records of macro evolution [19]. The non-stationary 
directional nature of the Tangled Nature model has been suggested to be an example of 
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generic properties of complex systems dynamics as observed in a number of very diverse 
phenomena by Anderson et al. [TO] . Here it was suggested that the directional gradual 
relaxation can be viewed as a release of a generalised strain originating in starting from 
a badly adjusted initial configuration. The Tangled Nature model suggests that in the 
case of ecosystems, selection and adaptation manage to direct the dynamics towards a 
selections of species better adjusted to coexist. 

In conclusion, it is useful to identify systems under the label "Complex Systems", 
since this indicates that the machinery of statistical mechanics can be applied to try 
to describe the system's dynamics and evolution. In ecology, there is an urgent need for 
"integrative" approaches that are non-stationary, and statistical mechanics can provide an 
initial mathematical framework, subject to modification and adaptation as one navigates 
deeper into the mysteries of complex systems. 
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